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ABSTRACT 
The  boundary  layer  for  a  swirling  flow  of  incompressible 
fluid  in  a  convergent  cone  has  been  investigated.   Theoretical 
calculation  with  boundary  conditions  more  appropriate  to 
physically  existent  situations  discount  the  existence  of 
"super  velocities"  within  the  boundary  layer.   Parallel 
experimental  investigations  demonstrate  an  interdependence  of 
core  and  boundary  layer  flows  which  precludes  the  maintenance 
of  initial  "ideal  flow"  conditions. 
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NOMENCLATURE 

u  velocity  in  x  direction  -  axial  velocity 

v  velocity  in  9  direction  -  tangential  velocity 

w  velocity  in  z  direction 

x  distance  along  the  wall 

z  distance  perpendicular  to  wall 

K  initial  ratio  of  tangential  to  axial  velocities 

U  free  stream  velocity  in  x  direction 

U1  derivative  with,  to  X  of  U 

V  free  stream  tangential  velocity 

V  derivative  with  respect  to  X  of  V 
a  acceleration  vector 

q  velocity  vector 

C  vorticity 

y  viscosity 

v  kinematic  viscosity 

p  density 

t  shear  at  the  wall 
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I.   INTRODUCTION 

The  determination  of  velocity  profiles  in  the  laminar 
boundary  layer  of  a  convergent  nozzle  has  been  approached 
analytically  in  several  different  manners.   Taylor  [1]  con- 
sidered the  case  of  a  predominant  tangential  flow  with  a 
secondary  axial  flow.   In  his  analysis,  Taylor  assumed  that 
the  tangential  flow  exceeded  the  axial  flow  to  the  extent 
that  the  axial  component   in  the  inviscid  core  could  be 
neglected.   After  applying  the  Pohlhausen  method  to  the 
mementum-integral  equations  and  obtaining  solutions  for  various 
axial  positions,  Taylor  concluded  that  a  condition  could 
arise  when  most  of  the  outflow  from  the  nozzle  was  fed  by  a 
boundary-layer  current  close  to  the  wall  of  the  swirl  chamber. 
The  radial  pressure  gradient,  which  held  the  particles  in  their 
circular  path,  acted  upon  the  retarded  boundary  layer  to  drive 
the  fluid  along  the  surface  to  the  exit. 

Wilks  [2]  developed  a  similar  approach  that  led  to  the 
possibility  of  a  'super-velocity'  existing  in  the  boundary 
layer  of  a  convergent  conical  nozzle.   This  'super-velocity' 
was  a  velocity  within  the  boundary  layer  that  exceeded  the 
free-stream  velocity.   Wilks  considered  the  flow  in  the  cone 
to  be  that  of  an  inviscid  core  with  a  thin  boundary  layer. 
The  inviscid  core  had  two  components  of  flow  -  a  uniform 
axial  flow  and  a  free   vortex  with  its  center  along  the  center 
line  of  the  nozzle. 

In  this  paper  an  analogous  approach  was  utilized.   How- 
ever, the  effect  of  the  boundary-layer  growth  upon  the  tangential 


and  axial  velocities  in  the  free  stream  flow  was  considered. 
Subsequently,  no  super-volocities  were  found  to  exist.   Further- 
more, experimental  results  prove  that  full  scale  effects 
predominate  in  this  problem  rather  than  just  boundary-layer 
effects. 
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II.   ANALYSES 

A.    WILKS  THEORY 

Using  the  Wilks  approach,  a  curvilinear  coordinate  system 

was  defined  with  the  x-axis  along  the  wall  and  the  z-axis 

perpendicular  to  the  wall  of  the  nozzle  (Figure  in  Appendix 

A) .   The  shape  of  the  nozzle  was  determined  by  the  surface  of 

revolution  of  an  arbitrary  curve  r(x).   Wilks  developed  an 

approximate  method  of  solution  using  the  momentum  equations 

(Appendix  A) 

3u     3u   v2dr      dU   V2dr      32u 

u  —  +w =U +v  — 5"  (1) 

3x     3z   r  dx      dx   r  dx      3z 

o 

3v      3v   uv  dr     3  v  /ox 

u  —  +w  —  + =v  — 5-  (2) 

3x      3z   r   dx     3z2 

|_  (r-u)  +  |_  (r-w)  =  0  (3) 

3x         3z 

and  the  corresponding  integrated  forms  of  the  momentum 
equations 

CO  00  .00 

i.  |_  (rfu(u-U)dz)  +  Hj  (u-U)dz  +  £-(  (V2-v2)dz      (4) 


=  -  v  (f|) 


o 


1   1_  (r2[  U(v-v)dz)  =  -  v  (5X)     ,  (5) 

r2  dx    Jo  3z  ° 


where  the  boundary  conditions  were 

2 

u  =  0;  -  v  ^  u  =  U^£  -  V2  £l;   on  z  =  0 
3z2    dx      r 

2 

u  =  U(x);  —  =  0;  — £  =  0;   when  z-h» 
3z       3z2 
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(6) 


and 


v  =  0;  d  v  =0;   on  z  =  0 
2 

(7) 


2 

v  =  V(x)  =  h-;    £l  =   0;  9  v  =  0;   when  z  -*•  « 

r(x)   8z      8z2 

The  distances  in  the  above  equations  were  non-dimensionalized 

by  defining  n  =  z  /  6  and  X  =  x  /  c,  where  6  was  the  local 

boundary-layer  thickness  and  c  was  the  distance  along  the 

generator  to  the  nozzle  to  its  apex. 

To  obtain  a  solution  of  the  momentum-integral  equations, 
it  was  necessary  to  assume  a  velocity  profile  within  the 
boundary  layer.   A  profile  that  would  also  satisfy  the  energy- 
integral  equations  should  produce  some  improvement  over  a 
solution  to  the  momentum- integral  equation  alone.   Wieghardt 
[3]  developed  such  a  method  with  the  use  of  two  free  parameters 
and  an  eleventh-degree  polynomial.   Applying  the  free  para- 
meters and  the  polynomial  to  the  u-velocity  profile  and  assum- 
ing that  the  v-velocity  profile  did  not  change  significantly 
with  X,  the  following  profiles  were  obtained: 

u/U  =  fx(n)  +  L1«f2(n)  +  L2*f3(n)  (8) 

where 

f  (n)  =  1  -  (1  -  n)8(l  +  8n  4-  36n2  +  120n) 

f2(n)  =  (1  -  n)8n(l  +  8n  +  36n2) 

fJn)  =  -(1  -  n)8n2(l  +  8n) 
J 
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and 


v/V  =  g(n)  =  2n  -  2n3  +  n4  (9) 


An  evaluation  of  the  axial  velocity  profile  at  the  wall  leads 
to  physical  connotations  for  the  shape  factors,  L-^  and  L^;  i.e 
Li  includes  the  wall  shear  stress  and  L2  involves  the  axial 
pressure  gradient. 


x06  _ -6  •  <5     dP 

yU  2    2c  yU  ^"  dX' 


Using  the  coefficients  Al,  A2,  A3  and  A4  (ratios  of  character- 
istic thicknesses)  as  defined  in  APPENDIX  A  and  introducing 
a  term 

G=0<£>=<#> 

the  momentum- integral  equations  became 

U«  fi'x    A2     r'^2 

__  (2  +  A1)  +   x  m   _  +  __y  A3          (10) 


—  + 

w2x 


rj_   UJ_     6^   A4 
r  +  u   +  62xy   G 


(11) 


As  a  method  of  simplifying  the  solution,  Wilks  assumed  that  a 
basically  linear  relationship  existed  between  the  momentum 
thickness  (<$2X)  and  the  boundary-layer  thickness  (6)  and 
between  the  mixed  momentum  thickness  ("^xv^  an<^  t^ie  koundary- 
layer  thickness.   The  above  equations  were  then  combined 
with  the  definition  of  G  and  its  derivative.   The  following 
relationships  were  obtained 
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r>  -  2.A4  -  (A2  -  A4>  -  r  r'  h  +    <v\2         A3   \  (12) 

G   ~  2  A4     (Al  -  1)      G  F~V  +  (U}    (1  +  Al)' 


A2  =  G  ~(1  +  Al)  +  A4  -  G '^1  "(S  2  A3  (13) 

U  r    U 

L  .  s  Jsi  -  Y*  i4  (14) 

z    2   lu     U^  r 

In  the  preceeding  analysis,  no  restriction  was  made  upon 
the  shape  of  the  convergent  nozzle.   The  flow  conditions  in 
a  cone,  as  well  as  many  other  shaped  nozzles,  do  not  meet  the 
requirements  for  a  similarity  solution  of  the  related  boundary 
layer  equations.   However,  in  these  cases  where  the  profiles 
vary  with  axial  position  in  a  non-linear  manner,  this  present 
approximate  method  of  solution  can  be  universally  used.   Wilks 
showed  that  a  solution  of  the  similarity  equation  could  be 
developed  for  any  one  point  along  the  surface  of  a  cone,  but 
that  the  coefficients  of  the  Falkner-Skan  equation  would  be 
different  for  each  point.   For  each  point,  he  noted  that  a 
close  agreement  existed  between  the  profiles  of  the  approximate 
method  of  solution  and  the  profiles  of  the  similarity  solution. 

For  the  particular  case  of  a  cone,  a  linear  relationship 
existed  between  the  radius  and  X.   The  axial  velocity  as  a 
functions  of  X  was  then  formulated  from  continuity,  and  the 
variations  in  the  tangential  velocity  were  derived  from  con- 
servation of  angular  momentum.   (The  details  of  these  functions 
are  given  in  APPENDIX  A.) 
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B.    NUMERICAL  SOLUTION  OF  WILKS  THEORY 

A  solution  of  the  momentum-integral  equations  would  lead 
to  the  values  of  the  shape  factors  and  thus  to  the  axial 
velocity  profile.   A  step-by-step  procedure  was  used  to  solve 
the  governing  equations  at  various  values  of  X.   Wilks 
considered  the  simplified  case  where  the  initial  boundary- 
layer  thickness  on  the  cone  was  zero.   With  this  assumption, 
G(0)  became  zero,  which  led  to  L?(0)  =  0.   Equation  13  then 
became  a  quadratic  in  L^(0).   The  solution  for  other  values 
of  X  continued  in  the  following  iterative  manner: 

(1)  Using  the  values  of  L-^(O)  and  L2(0),  the 
thickness  coefficients  were  calculated. 

(2)  The  velocities  were  then  evaluated. 

(3)  Using  equation  13,  L-.  was  calculated. 

(4)  Applying  equation  12,  the  value  of  G'  was 
determined. 

(5)  Using  the  first  two  terms  of  a  Taylor  series 
expansion,  G(X+dX)  =  G(X)  +  G' (X) »dX. 

(6)  L2  was  calculated  from  equation  14. 

(7)  The  values  of  the  thickness  coefficients  were 
recalculated  with  the  new  values  of  L-^  and  L2 
and  the  procedure  from  step  2  onward  was 
repeated  for  new  values  of  X. 

In  a  physical  system,  a  boundary  layer  could  easily  exist 
before  the  flow  entered  the  conical  region.   The  momentum- 
integral  equations  would  then  need  to  be  solved  for  a  non-zero 
initial  boundary-layer  thickness.   This  was  accomplished  by  a 
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slight  modification  of  the  Wilks  solution.   Using  the  solution 
of  the  quadratic  in  L-^(O)  for  zero  initial  boundary-layer 
thickness  as  a  predicted  value  of  L-^(O),  an  iterative  method 
would  correct  this  value  for  the  presence  of  an  initial 
boundary-layer  thickness.   The  method  pursued  was  the  following; 

(1)  An  initial  value  for  L-^(O)  was  assumed. 

(2)  The  velocities  were  calculated. 

(3)  G  and  thus  L2  were  determined  from  the  value 
of  the  initial  boundary-layer  thickness. 

(4)  Using  the  current  values  for  L  (0)  and  L  (0)  , 
the  thickness  coefficients  were  calculated. 

(5)  The  value  of  L-,  (0)  was  determined  according 
to  equation  13. 

(6)  The  calculated  value  of  L-,  was  compared  with 
the  assumed  value  of  Lt.   If  the  two  values 
were  not  sufficiently  close,  a  value  equal  to 
the  average  of  the  two  was  assumed  and  the 
procedure  from  step  4  onward  was  repeated. 

If  the  result  was  sufficiently  accurate,  L 
and  L~  were  calculated  for  other  values  of  X 
using  the  same  methods  as  for  a  zero  initial 
boundary-layer  thickness. 
Using  the  above  numerical  methods,  the  momentum-integral 
equations  were  solved  for  values  of  K  (swirl  factor)  of  one 
and  two.   The  resultant  profiles  were  compared  with  the  Wilks 
profiles  to  insure  the  correctness  of  the  numerical  method 
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employed.   The  solution  was  then  extended  for  values  of  K  up 
to  six  with  various  initial  thicknesses  and  initial  velocities. 
Several  profiles  that  were  obtained  in  this  manner  are  shown 
in  Figures  1-7.   The  solutions  indicated  that  the  degree  of 
super-velocity  increased  dramatically  with  K  and  X.   The 
value  of  the  momentum  thickness  became  negative  under  many 
circumstances,  and  for  high  degrees  of  super-velocity,  the 
displacement  thickness  also  became  negative.   Under  all  of 
these  conditions  the  boundary-layer  growth  continued  to 
increase  with  X.   The  numerical  analysis  deteriorated  for 
values  of  X  much  greater  than  0.3  in  most  of  these  solutions. 
There  was  also  a  failure  in  the  numerical  solution  for  ratios 
of  tangential  velocity  to  axial  velocity  in  the  range  of  three. 
It  was  believed  that  these  problems  were  caused  by  the  numer- 
ical method  employed  and  were  not  inherent  problems  of  the 
analytical  model.   The   quantity  (A2  -  A4)  could  be  very  close 
to  zero.   Since  one  of  the  terms  in  the  calculation  of  L, 
required  division  by  the  above  quantity,  overflow  in  the 
computer  could  easily  be  caused  when  this  quantity  approached 
zero.   A  different  final  arrangement  of  the  momentum-integral 
equations  would  avoid  this  problem  and  allow  for  a  solution 
at  higher  values  of  X. 
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FIGURE  1.   VARIATIONS  IN  PROFILE  WITH  K.   X=0 . 1   G(0)=0 
(a)  K=0  ,  (b)  K=l  ,  (c)  K=2  ,  (d)  K=3  ,  (e)  K=4 
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FIGURE  2.   VARIATIONS  IN  PROFILE  WITH  K.  X=0 . 1  G(0)=0.52 
(a)  K=l  ,  (b)  K=2  ,  (c)  K=3  ,  (d)  K=4 


19 


1.6 


1.4  - 


1.2  " 


1.0  - 


0.8  - 


0.6 


0.4 


0.2 


0.2 


0.4       0.6 


0.8 


1.0 


n 

FIGURE  3.   VARIATIONS  IN  PROFILE  WITH  K.  X=0.2  G(0)=0.52 
(a)  K=l  ,  (b)  K=2  ,  (c)  K=4 
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FIGURE  4.   VARIATIONS  IN  PROFILE  WITH  X.  K=l  G(0)=0.52 
(a)  X=0.00  ,  (b)  X=0.05  ,  (c)  X=0.10,  (d)  X=0.15 
(e)  X=0.20  ,  (f)  X=0.25 
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FIGURE  5.   VARIATIONS  IN  PROFILE  WITH  X.  K=2  G(0)=0.52 
(a)  X=0.00  ,  (b)  X=0.05  ,  (c)  X=0.10  ,  (d)  X=0.15 
(e)  X=0.20  ,  (f)  X=0.25 
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FIGURE  6.   VARIATIONS  IN  PROFILE  WITH  X.  K=4  G(0)=0.52 

(a)  X=.00  ,  (b)  X=.05  ,  (c)  X=.10  ,  (d)  X=.15  ,  (e)  X=.20 
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FIGURE  7.   VARIATIONS  IN  PROFILE  WITH  G.  K=2  X=0.2 
(a)  6(0)=0.00  ,  (b)  6(0)=0.10  ,  (c)  6(0)=0.15 
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C.    PRESENT  THEORY 

The  Wilks  approach  to  the  momentum  equation  assumed  that 
the  tangential  velocity  at  the  outer  edge  of  the  boundary 
layer  could  be  described  as  A/r.   This  was  only  appropriate 
if  the  boundary-layer  thickness  was  very  small  compared  to 
the  value  of  r.   The  numerical  solution  of  the  Wilks  method 
led  to  a  boundary-layer  growth  that  could  be  of  the  same 
order  as  r  for  flow  in  air  where  r  was  about  one  inch  and  the 
initial  axial  velocity  was  one  foot  per  second.   In  these 
cases,  the  radius  minus  the  boundary-layer  thickness  should 
be  used  to  determine  the  tangential  velocity  at  the  outer 
edge  of  the  axial  boundary  layer  rather   than  using  the  radius 
alone.   Thus  V  =  A/(r  -  6). 

The  introduction  of  this  new  value  of  V  into  the  boundary 
conditions  of  the  momentum  equations  led  to  a  different  form 
of  the  v  momentum-integral  equation.   The  derivation  of  the 
equation  is  shown  in  APPENDIX  B.   This  approach  required  the 
definition  of  a  new  thickness  coefficient  A5  (APPENDIX  B) . 
The  v  momentum-integral  equation  then  became 


Jr'  ,  V       A5    [6*         r'\l    A4 


*L  +  LL_   §-'     ■  Jr'  ■  V'  ■     A5    '6'    r'H  _  A4    (15) 
r    U    6 


In  the  above  equation,  Wilks  assumption  of  a  basically  linear 
relationship  between  the  characteristic  thicknesses  was  used. 
The  terms  in  the  braces  are  added  to  the  Wilks  v  momentum- 
integral  equation  by  the  new  assumption,  vix.,  V  =  A/(r  -  6). 

If  the  boundary  layer  is  thick  enough  to  effect  the 
boundary  condition  for  the  tangential  flow,  the  calculation 
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of  the  axial  velocity  should  also  take  this  boundary-layer 
growth  into  consideration.   In  using  the  conditions  of  conti- 
nuity of  flow,  the  radius  minus  the  displacement  thickness 
should  be  used  rather  than  the  radius  alone .   The  displacement 
thickness,  by  its  definition,  is  the  distance  the  potential 
flow  has  been  displaced  by  the  presence  of  the  boundary  layer. 
The  effect  of  this  consideration  would  be  to  increase  the 
axial  core  velocity  as  the  displacement  thickness  grew. 

D.    NUMERICAL  SOLUTION  OF  PRESENT  THEORY 

In  order  to  calculate  U  and  its  derivative,  U' ,  the  value 
of  the  displacement  thickness  and  its  derivative  were  needed. 
The  ratio  of  the  displacement  thickness  to  the  boundary-layer 
thickness  was  one  of  the  integrals  used  in  the  calculation 
of  the  thickness  coefficients.   The  boundary-layer  thickness 
at  any  point  could  easily  be  calculated  from  the  value  of 
G  once  the  kinematic  viscosity  and  length  of  the  nozzle  to 
its  apex  had  been  specified.   The  calculation  of  the  derivative 
of  the  displacement  thickness  was  not  as  direct.   Wilks  based 
his  solution  on  the  assumption  that  a  basically  linear  relation- 
ship existed  between  the  momentum  thickness  and  the  boundary- 
layer  thickness.   An  extension  of  this  assumption  would  be 
to  presume  a  basically  linear  relationship  also  existed  between 
the  displacement  thickness  and  the  boundary-layer  thickness. 
Thus  the  ratio  of  the  displacement  thickness  to  its   derivative 
should  be  equal  to  the  ratio  of  the  boundary-layer  thickness 
to  its  derivative. 
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With  the  new  equation  for  the  v  momentum  integral,  the 
numerical  solution  proceeded  in  the  following  manner: 

(1)  Using  the  initial  boundary-layer  thickness, 
the  values  of  G(0)  and  L2  (0)  were  calculated. 

(2)  A  predicted  value  for  L, (0)  was  assumed  from 
the  Wilks  solution  for  zero  initial  boundary 
layer  thickness. 

(3)  The  distance  along  a  flat  plate  that  would  be 
required  to  produce  a  boundary-layer  thickness 
of  the  same  size  as  the  initial  boundary-layer 
thickness  was  calculated.   Using  this  distance, 
the  derivative  of  the  boundary-layer  thickness 
on  a  flat  plate  at  this  point  was  calculated. 
From  this  estimate  of  the  derivative  and 
assuming  an  initial  value  for  Dl  from  the  Wilks 
solution,  the  initial  value  of  the  derivative 
of  U  was  evaluated. 

(4)  Using  the  current  values  of  L]_  and  L2  /  the 
values  of  the  thickness  coefficients  were 
calculated. 

(5)  The  value  of  6'/<5  was  calculated,  where 

6'    ,A4    U'    r'  (0  A5       V\  . 

T~   "  lG    U     F~  U  -  r/6  -  1'         V  )/ 

(1  +'  ^7^'  (16) 
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(6)   The  new  value  of  L,  was  calculated  from 


Lx  =  D2*G« 


—  (2  +  Al)  +  £!  (1  -  3L  A3)  + 
U  r       U2 


H 


(17) 

(7)  The  calculated  value  of  L-^(O)  was  compared  with 
the  assumed  value.   If  the  two  values  were  not 
sufficiently  close,  the  thickness  coefficients 
and  velocities  were  re-evaluated  until  the 
values  of  L, (0)  for  two  consecutive  iterations 
were  sufficiently  close. 

(8)  The  value  of  G'  was  determined  from 

(9)  The  values  of  G  and  L?  for  the  next  point  were 
then  calculated. 

(10)  The  step-by-step  evaluation  of  the  momentum- 
integral  equations  was  continued  for  new  values 
of  X  by  returning  to  step  4  and  remembering 
that  step  7  is  used  only  when  X=0 . 

U  and  V  as  functions  of  X  are  given  in  APPENDIX  B. 

Using  this  numerical  technique,  the  momentum-integral 
equations  were  solved  for  various  values  of  initial  boundary- 
layer  thicknesses,  initial  velocities,  and  K.   Since  this 
method  required  the  division  by  G  in  the  calculation  of  equa- 
tion 16,  the  analysis  was  not  valid  for  the  case  where  the 
initial  boundary-layer  thickness  was  zero.   However,  for  very 
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small  values  of  G  other  than  zero,  there  was  no  restriction 
on  the  initial  boundary-layer  thickness.   The  solution  of 
these  equations  showed  no  sign  of  super-velocity.   In  several 
cases  the  step-by-step  evaluation  continued  smoothly  for  values 
of  X  up  to  0.6.   The  boundary-layer  growth  reached  a  maximum 
value  when  X  was  in  the  range  of  0.25  and  then  decreased  with 
X.   This  numerical  analysis  appeared  to  lose  its  validity  when 
the  momentum  thickness  became  negative,  which  occurred  at  a 
much  higher  value  of  X  then  in  the  Wilks  solution.   At  this 
value  of  X  where  the  momentum  thickness  became  negative,  dis- 
continuities appeared  in  the  boundary-layer  growth.   Figures  8-12 
show  a  comparison  between  the  profiles  of  the  Wilks  solution, 
as  modified  for  initial  boundary-layer  thickness,  and  the 
profiles  obtained  by  the  above  method. 

E.    DISCUSSION  OF  THEORETICAL  RESULTS 

Wilks  based  his  solution  of  the  flow  in  the  boundary  layer 
of  a  convergent  nozzle  on  several  assumptions.   He  first  con- 
sidered the  inviscid  core  as  a  free  vortex  with  uniform  axial 
flox.   It  was  then  assumed  that  the  boundary  layer  would  be 
thin  enough  to  prevent  any  effect  of  its  growth  upon  the  core 
velocities.   Thus  the  core  velocities  were  dependent  only  upon 
the  geometric  configuration  of  the  convergent  nozzle.   Wilks 
also  assumed  that  a  linear  relationship  existed  between  the 
various  characteristic  thicknesses  of  the  boundary  layer.   With 
these  assumptions,  the  solution  of  the  momentum-integral 
equations  led  to  the  existence  of  a  super-velocity  in  the 

boundary  layer. 
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For  low  velocity  air  flow  (required  to  maintain  the 
laminar  conditions)  the  boundary-layer  growth  was  not  negli- 
gible.  By  including  this  growth  in  the  calculations  of  the 
boundary  conditions,  the  solution  no  longer  led  to  super- 
velocities.   Increased  swirl  only  decreased  the  ratio  of  the 
displacement  thickness  to  the  boundary-layer  thickness 
slightly. 
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III.   EXPERIMENTATION 

A.    DESIGN  OF  APPARATUS 

The  analytical  approach  assumed  that  the  flow  entering 
the  convergent  nozzle  was  a  combination  of  a  free  line-vortex 
and  uniform  axial  flow.   To  test  the  theory,  an  apparatus  was 
required  to  approach  this  assumed  velocity  distribution.   There 
are  several  methods  of  introducing  a  fluid  into  a  chamber  and 
causing  a  swirling  motion.   However,  this  tangential  velocity 
would  have  to  approximate  the  conditions  of  a  free  vortex,  that 
is  the  tangential  velocity  should  be  inversely  proportional  to 
the   radius.   Several  methods  were  considered,  including  rotating 
a  section  of  the  swirl-producing  chamber,  drawing  the  fluid  into 
the  chamber  from  a  stagnation  condition  while  using  vanes  to 
create  the  tangential  motion,  and  rotating  blades  inside  the 
chamber.   The  approach  chosen  was  to  introduce  the  fluid  through 
one  tangential  port  with  variable  cross-section.   In  tests  of  a 
Hilsch  tube,  Lay  [4]  used  a  single  circular  port  to  introduce 
a  jet  into  the  chamber.   This  was  the  only  source  of  fluid  in 
his  tests,  and  he  reported  the  velocity  profiles  that  were 
produced  in  this  manner.   However,  he  did  not  indicate  the 
radial  profile  of  the  velocities  in  a  region  close  to  the  wall. 
The  results  indicated  that  the  tangential  flow  might  closely 
approximate  the  conditions  for  a  free  vortex  over  most  of  the 
radius,  departing  from  these  conditions  in  the  center  portion. 
The  distribution  of  the  axial  velocity  in  Lay's  tests  suggested 
that  the  superposition  of  pipe  flow  upon  this  flow  could  result 
in  a  nearly  uniform  axial  flow. 
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The  source  of  fluid  for  this  apparatus  was  a  compressed 
air  tank.   To  insure  laminar  pipe  flow  of  the  air  at  speeds 
of  at  least  two  feet  per  second,  a  pipe  diameter  of  two  inches 
was  chosen.   Axial  pipe-flow  could  be  produced  prior  to  entry 
of  the  axial  flow  into  the  swirl  chamber.   The  shape  of  this 
axial  flow  could  be  varied  by  adjusting  the  axial  position  of 
a  honeycomb  f low-straightener .   If  profiles  different  from 
pipe  flow  were  required  to  achieve  the  uniform  axial  profile 
at  the  entrance  of  the  cone,  the  flow  could  be  adjusted  by 
only  passing  air  through  selected  annular  portions  of  the 
honeycomb  f low-straightener. 

The  tangential  flow  was  introduced  through  a  rectangular 
slot  in  the  wall  of  the  swirl  generator.   One  wall  of  this 
slot  was  tangent  to  the  inner  surface  of  the  chamber.   The  slot 
was  0.25  inches  high  and  the  width  could  be  varied  from  0.3 
to  1.25  inches.   This  variable  cross-section  permitted  an 
adjustment  in  the  ratio  of  the  velocity  to  the  quantity  of  mass 
introduced  by  the  tangential  port.   This  flexibility  in  the 
system  was  added  to  insure  that  the  proper  profiles  and  values 
of  K  could  be  obtained.   Tests  were  then  performed  to  determine 
the  axial  distance  required  for  the  development  of  the  appro- 
priate flow  conditions.   These  preliminary  tests  indicated  that 
a  range  between  one  and  two  feet  was  required  for  the  development 
of  the  tangential  jet  into  a  free  vortex. 

Lay's  tests  indicated  that  there  would  be  a  region  close 
to  the  wall  where  the  tangential  flow  would  be  retarded  by 
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viscosity.   An  axial  combination  of  pipe  flow  and  the  tangen- 
tially  introduced  flow  would  also  have  a  region  close  to  the 
wall  where  the  conditions  for  uniform  axial  flow  would  not  be 
met  due  to  the  viscous  forces  present.   With  these  two  boundary- 
layer  growths  under  consideration,  the  outer  annular  portion 
of  the  flow  was  bled  off,  with  only  the  "flat"  inner  portion 
of  the  axial  flow  entering  the  cone.   This  bleed-off  was 
accomplished  by  placing  a  cone  of  initial  radius  0.77  inches 
at  the  end  of  the  tube  which  had  an  inside  radius  of  one 
inch.   The  position  of  the  cone  relative  to  the  end  of  the 
tube  could  be  varied  in  order  to  alter  the  amount  of  flow  bled 
off  into  the  atmosphere.   The  entry  region  into  the  cone 
consisted  of  a  one-inch  straight  section.   The  purpose  of  this 
section  was  to  insure  the  continuation  of  the  proper  flow 
characteristics  after  bleed-off.   (See  Figure  25) 

The  cone,  a  piece  of  machined  plexiglass,  featured  a 
semi-angle  of  4.4  degrees.   The  distance  along  the  generator 
of  the  cone  to  the  apex,  termed  c  in  the  analysis,  was  ten 
inches.   The  truncated  cone  used  in  the  test  apparatus  extended 
over  0.35  c.   Thus,  tests  could  easily  be  performed  over  the 
region  where  X  varied  from  0.0  to  0.3.   Two  hot-wire  probes 
were  introduced  from  the  exhaust  end  of  the  truncated  cone. 
One  configuration  was  sensitive  to  axial  and  radial  components 
alone.   The  other  probe  could  detect  tangential  and  radial 
velocities.   It  was  assumed  that  the  radial  component  of  flow 
was  negligible,  especially  close  to  the  wall.   The  traversing 
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mechanism  moved  in  two  directions,  allowing  for  motions  in 
both  the  x-direction  and  the  z  direction. 

B.    EXPERIMENTAL  RESULTS 

The  development  of  flow  conditions  that  would  approximate 
the  desired  conditions  was  basically  a  trial-and-error  pro- 
cedure.  By  adjusting  the  slot  width,  the  development  length 
for  the  tangential  flow  (by  changing  the  length  of  pipe  between 
the  swirl  generator  and  the  nozzle),  the  bleed-off  level,  and 
the  radial  distribution  of  the  axially  introduced  flow,  it  was 
possible  to  achieve  several  flow  conditions  that  approximated 
the  analytical  boundary  conditions  at  the  entrance  of  the  cone. 
However,  in  all  the  cases  tested,  the  velocity  distribution  in 
the  inviscid  core  changed  drastically  with  X.   The  radial  dis- 
tribution of  the  axial  and  tangential  velocities  for  several 
values  of  X  with  swirl  ratios  (K)  of  0.6,  0.8,  and  0.95  are 
shown  in  Figures  13-24.   These  profiles  indicate  several 
different  flow  relationships.   For  the  case  where  K=0 . 6 ,  the 
degree  of  swirl  had  very  little  effect  upon  the  axial  flow  as 
it  moved  through  the  cone.   The  uniform  axial  core  was  gradually 
replaced  by  a  flow  that  resembled  the  distribution  in  pipe  flow, 
with  most  of  the  flow  occuring  in  the  central  portion,  and  the 
perimeter  retarded  by  viscous  forces.   With  the  slight  increase 
to  K=0.8,  there  was  a  great  change  in  the  radial  distribution 
of  the  axial  flow.  In  this  case,  the  free  vortex  seemed  to  be 
maintained  for  all  test  stations.   At  X=0.3,  this  tangential 
motion  had  thrown  a  great  deal  of  the  fluid  into  the  peripheral 
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portion  of  the  cone,  leading  to  a  high  axial  velocity  near  the 
wall.   However,  this  high  velocity  was  still  in  the  core 
(z=0.2  inches)  and  not  in  the  boundary  layer.   In  the  tests 
where  K=0.95/  the  distance  over  which  the  tangential  flow 
appeared  as  a  free  vortex  was  less  then  in  the  previous  cases. 
This  smaller  distance  resulted  in  several  changes  in  the 
tangential  velocity  with  X.   As  a  result  of  this  changing 
tangential  profile,  there  were  several  in  the  axial  flow.   In 
sections  where  the  tangential  flow  appeared  as  a  forced  vortex, 
there  was  one  form  of  axial  flow,  where  high  speed  occurred 
away  from  the  wall,  while  in  annular  sections  where  the  tan- 
gential flow  resembled  a  free  vortex,  the  high  speed  axial  flow 
occurred  nearer  the  wall. 

C.    DISCUSSION  OF  EXPERIMENTAL  RESULTS 

In  both  of  the  analytical  approaches,  the  inviscid  core 
was  considered  a  combination  of  uniform  axial  flow  and  a  free 
vortex.   In  the  test  cases,  the  axial  flow  at  the  entrance  of 
the  nozzle  was  within  five  per  cent  of  uniform  axial  flow.   The 
tangential  flow  approximated  a  free  vortex  over  most  of  the 
radius  in  the  entrance  region.   However,  this  flow  could  not 
be  maintained  down  the  cone.   Swirl  induced  a  high  velocity 
axial  flow  near  the  wall  when  the  tangential  flow  was  in  the 
form  of  a  free  vortex.   This  region  of  high  axial-velocity  flow 
was  far  enough  away  from  the  wall  (a  distance  greater  than  0.15 
inches)  to  call  it  part  of  the  inviscid  core.   Thus  swirl  does 
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indeed  cause  changes  in  the  axial  flow,  but  these  changes  occur 
in  the  core  and  not  in  the  boundary  layer  as  Wilks  suggested. 
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IV.   CONCLUSIONS  AND  RECOMMENDATIONS 
Previous  work  in  this  area  considered  the  effect  of 
swirling  flow  upon  the  velocity  profiles  in  the  boundary 
layer.   In  this  paper,  the  effect  of  the  boundary  layer  growth 
upon  the  core  velocities  was  utilized  in  the  analysis.   How- 
ever, the  experimental  tests  showed  that  this  consideration 
was  not  enough.    The  swirl  also  affected  the  axial  core  flow, 
and  the  resultant  changes  are  much  more  drastic  than  any  varia- 
tion in  the  boundary  layer  flow.   Rather  than  applying  the 
Navier-Stokes  equations  to  the  boundary  layer  alone,  it  would 
be  much  more  appropriate  to  apply  them  to  the  whole  flow  field. 
In  addition,  the  developed  theory  for  the  boundary  layer  assumed 
that  a  basically  linear  relationship  existed  between  the 
various  characteristic  thicknesses.   Since  these  thicknesses 
are  functions  of  the  shape  factors,  which  are  in  turn  functions 
of  X,  it  would  seem  mandatory  to  extend  this  analysis  without 
these  assumptions  of  linear  relationships.   However,  in  order  to 
complete  such  an  analysis,  more  information  needs  to  be  known 
about  the  effect  of  swirl  upon  the  core  velocities. 
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APPENDIX  A 
DERIVATION  OF  EQUATIONS  ACCORDING  TO  WILKS 
The  boundary  of  the  convergent  nozzle  is  formed  by  sur- 
face of  revoltuion  of  a  curve  r(x),  where  x  is  a  distance 
measured  along  the  wall.   The  distance  from  the  wall  is  z. 
The  potential  flow  has  the  two  components,  V  and  U  as  shown 
below 


Acceleration  in  the  orthogonal  curvilinear  coordinates  is 
defined  as 


» -  i  - 1 + '  w*2>  -  5  x  ? 


In  the  i .  direction  this  becomes 

1 


'dq±         q±    dqL         q2    dq±         q3    8qi 

Tt~  +  H7   9u7  +  h~   au~  +  h~  9u~ 
11  2        2  3        3 


h„     \dun)     +        h_      \8u    / 


2  2 

h2      ^dui' 
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3 
9h 


ho      Vsu^/i 


If   u^  =  x   /  .  1*2  -  6   /   U3  =  z 

<2l  "  u  f      <22  =  v  1      S3  =  w 
h^  -  1   /   h2  =  r   ,   I13  =  1 

the  acceleration  can  be  written 

8u   .  •  3u ,  V  3u      3u   :  V2  3r 

x   at     3x   r  36      3z    r   ax 

3w     8w   V  aw      aw    V2_   3r 

az  =  at  +  u  3x  +  r  ae  +  w  az  "  r  az 

3v     8v   v  av      av   vu  3r 
ae  "  at  +  u  8x  +  r  36  +  w  3¥  +  ~F~  3x. 

For  axisymmetric  flow,  all  of  the  partial  derivatives  with 
respect  to  the  angular  displacement  are  zero.   Since  velocities 
in  the  direction  perpendicular  to  the  wall  (the  z  direction) 
may  be  considered  small  compared  to  the  other  velocities, 
the  momentum  equation  in  the  z  direction  may  be  neglected.   By 
an  order  of  magnitude  analysis,  the  boundary-layer  equations 
for  steady  laminar  flow  with  constant  properties  become 


2  ,      ,  _     _  .2 

L 
2 


■  --lE_y^&r__13P  +  v2  £fH 


ax      az    r  dx     p  8x      3z 


3v      3v   uv  dr     1  3P     2  32v 
u  —  +  w  —  + =  -  —  —  +  v 


ax      az    r  dx     p  30      3z 


2 


I;  (r-u)    +     I-  (r-w)   =  0 

with  the  boundary  conditions 
u  =  v=0      onz  =  0 
u  *  U(x)   and  v  -»■  V(x)  =  A/r(x)     as  z  -*■  » 
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With  the  above  boundary  conditions,  the  pressure  gradient  in 
the  x-direction  can  then  be  written 

-  I   i£  .   u  ^  -  V2 '  £l 
p  3x      dx       r 

Integrating  the  u-momentum  equation  across  the  boundary  layer 
gives 


f  f      3u     8u   vz 


dr     dU 

— —  -  IT      4- 

dx    u  dx  + 


v2|l)dz=-v(|H)o 


from  continuity 

f'2 


W  =  - 


o  3x 


(r*u)  dz 


Thus 


J0\   9x   r  8z 


4-  (r-u)  dz  -  v^  £_ 
0  9x  r 

v2  |l)az  =  -  v(|f) 


-   u 


dU 
dx 


The  second  term 


can  be  integrated  by  parts. 


m  = 


=   f0  lx  (r'U)  d: 


dp  = 


3i"  dz 


dm  =  -7T—    (r*u)  dz 

oX 


p   =  u 
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oo      ,00 


m  p 


-J 

*  0 


p  dm  =  u 


■5—  (r»u)dz 
9x 
0 


u  j^    (r«u)dz 


00  .00 


=  U    L_  (r»u)dz   - 
8x 


u  _  (r-u)dz 
8x 


The  u  momentum-integral  equation  then  becomes 

f  fu  ^H  -  H  i_  (r-u)  +  £  <L_  (r-u)  -  U  ^i\   d: 
J0V   9x   r  3x         r  3x  dx  J 


(V 


2  _  „2.  r 


vz)  ±_  dz  = 

r 


-  v 


(It). 


The  first  integral  can  be  written 


f  (u  |»  -  UZH.  |r  |«  +  u ■  |ri.  0  du\  dz 

/  \   9x    r   *   9x     3x/     dx J 

f°Y?  «  9u   n  9u    (U-u) -u  dr    „  du\  , 

-    z  u  -s —  -  u  -r —  -  - - -5—  -  u  3 —  J  d: 

J  0\  dx  dx  r    dx     dx/ 


Since 


•  s 


9uU 
9x 


u 


au 

dx 


the  integral  becomes 


(°°f 
J  cA 


5u 
9x 


2 


9(uU) 
dx 


(u-U)    *E  +    ^U)'u  **)     d: 
dx  r  dx  J 


b  1L    fr  u-  (u-U)    dz  ] 

r    9x     V       Jo  / 


dU 
dx 


(u-U)    dz 


The   integral    form  of   the  v-momentum  equation   is 


f  (u  |v  +  „ ■  |v  +  vvy  dr  \ 
J    I       9x  9z  r      dx  J 


dz   =   -   v 


»). 
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Introducing  the  continuity  equation,  this  becomes 


0  a  -  i 


l0  3x 


u*v  dr  \ 
r  dx } 


(f-u)  dz  +  HIV  §£  dZ  =  -  v 


(£). 


Performing  an  integration  of  the  second  term  by  parts 
•z 

o 


m 


-i: 


fL.  (r*u)  dz 

8x 


dp   =  I  IX  dz 

r  3z 


dm  =  7t—  (r«u)  dz 

9x 


_  v 


m  p 


■^  (r*u)  dz  -    p^  (r-u)  d: 
o  '  o 


The  left  hand  side  of  the  momentum-integral  equation  may  now 
be  written  as 


-CO   - 


|_{r.u)  +^(r.u,  *!£2g, 


) 


=  frf9(v'u)  -  v  23^  -  v  2S  +  v  1H  -  (v-v)  u 

j0Lv8x      ~   8xJ     3x     3x     r 

+  Xm  dr]  dz 
r  dxj 

=  f"f|_  (vu)  -  V  |H  -  Vm.  |r  +  2Iu^  dr  \  &z 
j  I  8x  dx  r      dx  r        dx J 


8r 
8x 


Since 


the  equation  can  be  written 

(  (<L_  (vu  -  Vu)  +  llH  (v-V)£*E^dz  =  -  vflX\ 
J0V9x  r       dxj         Uzj 
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Thus  the  v  momentum-integral  equation  is 


1_  d_ 
r2  dx 


u  (v-V)  dz   =  -  v 


m\ 


with  the  boundary  conditions 
v  =   0 


on   z  =  0 


1_Z  =   0 
8z2 


V        = 

V(x)    = 

A 
r(x) 

■^ 

8v 

dz      ~ 

0 

~2 
8    v 

J 

when  z  -♦  °° 


8z' 


=   0 


The  u  momentum-integral  equation  is 


1  d_ 
r  dx 


f   u  (u- 


U)  dz    + 


dU 
dx 


(u-U)  dz 


+  1  dr 
r  dx 


r  <v2-v2)  az = - v  in) 


with  the  boundary  conditions 
u  =   0 


-  v 


92u 

8^2 


=   U 


dU 


V2  dr 


dx   r  dx 


►  on   z  =  0 
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.  u   =   U(x) 

|H  =  o 

dz  u 


92u 


►  when  z  -»■ 


~2   =   0 
3z 

In  order  to  obtain  a  solution  to  the  moment urn- integral 
equations  using  the  approximate  method  of  von  Karman,  it  is 
necessary  to  assume  suitable  velocity  profiles.   A  direct 
analogy  with  the  Pohlhausen  method  would  lead  to 

u/U   =   f(n,Lx)      v/V  =   g(n,L2) 
where  n   =   z/6 (x) .   The  continuity  equation  leads  to  a  stream 
function 

$  =  rsU»h(x,n) 
The  axial  velocity  can  then  be  written 

u  =  U  A-  h(x,n) 
dn 

If  the  tangential  velocity  is  of  the  form  v  =  V*g(x,n)  , 

the  momentum  equations  may  be  written 

h,,r  +  a-h-h"  +  b'(l-h'2)  +  d-  (1-g2)  -  0 
g"  +  a«h»g'   =  0 
where  a,  b,  and  d  are  constants  as  defined  by  Wilks.   Boundary 
conditions  require  that  g"  evaluate'  at  the  wall  is  zero.   A 
differentiated  form  of  the  v-momentum  equation  gives 

g"'  +  a» (h' «g'  +  h«g")   =   0 
Since  both  hl  and  g"  evaluated  at  the  wall  are  zero,  it  is 
apparent  the  g" '  evaluated  at  n  =  0  must  also  be  zero.   This 
would  require  the  second  and  third  derivatives  of  g  to  be 
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independent  of  x.   With  this  understanding,  the  two  parameters 
can  be  used  to  obtain  the  axial  velocity  profile  while  the 
tangential  profile  is  considered  dependent  upon  n  alone. 
Thus 

u/U  =  f  (n,L  ,L  )   and   v/V  =  g(n) 
Introducing  a  dimensionless  length  X=x/c  (c  being  the  length 
of  the  cone  along  the  generator) ,  the  u  momentum-integral 
equation  becomes 


^Mr,u2-6io 


f  (f-l)dn 


1  dU 
+  c  dX 


U-6  [ 

J  n 


(f-l)dn 


+  _i_  V2  <£.-& 
c-r    3x 


(l-g2)dn  =  -  v  U  f • (0) 

o 


with  the  boundary  conditions 


f 

= 

0 

«2 

= 

U   dU   . 
c   dX 

V2      dr 

cr  dX 

on  n  =  0 


f   =   1 

f '   =   0 
f" '  =   0 

The  v  momentum-integral  becomes 


on  n  =  1 


—. -,   —   r2-U-V-6     f(l-g)dn   =  -  V-g1 
•r   dX  \  J0         /     6 


with  the  boundary  conditions 


=   0 
=   0 


*     on  n  =  0 


(0) 
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g    =   1 

g'   =0        ^  on   n  =  1 

g"   =   0 

Here  it  was  assumed  that  the  boundary  conditions  could  be 
transferred  from  infinity  to  the  boundary-layer  thickness, 
that  is  where  n=l,  without  appreciable  error. 

In  an  attempt  to  obtain  an  improvement  in  the  basic 
approximate  method  of  solution,  a  profile  can  be  obtained  that 
will  also  satisfy  the  requirements  of  the  energy-integral 
equation.   Weighardt  [3]  developed  such  a  method  using  two- 
parameter  profiles.   This  profile  can  be  represented  with  the 
aid  of  an  eleventh  order  polynomial.   Thus 

u/U  =   fx(n)   +  Lx-f2(n)   +  Iyf^n) 

where 

f  (n)   =   1  -  (l-n)8(l  +  8n  +  36n2  +  120n3) 

f2(n)   =   (l-n)8n  (1  +  8n  +  36n2) 

f3(n)   =   -(l-n)8n2(l  +  8n) 

Evaluating  the  functions  at  the  wall  (n=0)  yields 


/9u\  _   yU  (df\        _      uU  , 


thus 

Li  -  w 
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Also 


VII 


U  dU   V    dr 


-  "7-  f"  = —   at   n=0 

6'"      c  dX   cr  dX 


where  f" (0)  =  -2  L2 


This  leads  to 


L0  = 
2    v 


U6' 


L.    f  Hi  [V\2   rj 
2c   U  *\v)        r  J 


The  integrals  in  the  momentum  integral  equations  can  be 
expressed  as  various  characteristic  thicknesses. 


6n     =6 
lx 


62x   =  « 


60    =6 
2xy 


o 
1 

0 

'1 

.  0 


(l-f)dn 


f  (l-f)dn 


displacement  thickness  in  the 
x  direction 


momentum  thickness  in  the  x 
direction 


(l-g)dn      mixed  momentum  thickness 


With  the  definition  of  a  new  thickness 

.1 
o 
the  momentum-integral  equations  can  now  be  written 


6     (l-g2)dn 
J  o 


_L_d   (r-U2-60  )  +I^U  6lv-  JL  V2  gr  6   =  v  U  f,  (Q) 
cr  dX  u      2x'    c  dX    lx  c«r    dX   4     % 


_i_d_  (r2.u.v.a2xy)  -»  V.g'(O) 
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Differentiating  the  appropriate  terms,  these  equations 
simplify  to 

r1   U'        ,    6SX   vc        r'  V2 

—  +  —  (2  +  Al)  +  -**  =  2  A2  + 7  A3 

r    U  62x   u»6       r  U 

r'    U«     6'      vc 

—  +  —  +  — =£*  =  A4 

r    U     62xy    U-62 


where 


i: 


(l-f)dn 

6lx    D1 

AT  =   — t =    -*-x   = 

I1  6?v     D2 

f(l-f)dn      2x 


A2  =   ■  f(Q)    J  *1 

'1  D2 

f  (l-f)dn 


i: 
i: 


(l-g2)dn 


A3  =   _2 =  PI 

D2 


C 


f  (l-f)dn 


g"(Q)    2 

f(l-g)dn 


An  integration  of  the  polynomials  in  the  above  integrals  led 
to 

Dl   =   0.3333333  -  0.0454545  L1  +  0.006060  L 
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D2   =   0.0740342  +  0.0149190  L±    -  0.00118583  L 

-  0.0045481  L  2  -0.00008286  L22  +  0.0012013  L1L2 

D3   =   0.08461538  +  0.02435846  L,  -  0.00293040  L 

1  2 

D4   =   0.4174603 
Wilks  made  the  assumption  that  a  basically  linear  relationship 
existed  between  the  boundary-layer  thickness  and  the  momentum 
thickness.   This  assumes  that  D2  does  not  vary  significantly 
with  X.   A  basically  linear  relationship  between  the  mixed- 
momentum  thickness  and  the  boundary-layer  thickness  was  also 
assumed. 

With  the  introduction  of  a  dimensionless  parameter 

the  foregoing  equations  can  now  be  written 


§1  +  £1  (2  +  Al)  +  «1  =  >£   +  El  (^)2  A3 


r    U   v      '         6    G    r 


r  +  U  +  6   ~  G 


Now 


£1       iL    yj_ 

G   "  2  6   +  U 


Combining  these  three  equations  yields 


2 

A2  =  G  El  (1  +  Al)  +  A4  -  G  El(Y) 
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In  addition  to  the  above  relationships,  the  previously  derived 
equation  for  the  shape  factor  L?  is  also  important,  viz., 

2  _, 


L2  =  f(FGe  f] 


The  following  characteristics  must  be  considered  for  the 

perfectly  conical  nozzle. 

r(X)  =  r  (1  -  X) 
o 

By  continuity 

U  -r2'7r   =   U(X)  r2»(l  -  X)2»tt 
o  o  o 

•  u^ 

U(X)  = 


(1  -  X)2 


Conservation  of  angular  momentum  requires  that 

V  «r   =  V(X) -r  • (1  -  X)  =  A 
o  o         o 

If  V   =  K«U 
o      o 

v(x)  -  K'uo 


(1  -  X) 

Furthermore , 


Uj_       2 
U   '    (1  -  X) 


and 

r' 


(l  -  x) 
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APPENDIX  B 
DERIVATION  OF  PRESENT  THEORY 

The  concept  of  a  free  vortex  in  the  inviscid  core  requires 
that  the  tangential  velocity  be  inversely  proportional  to  the 
radius.   In  the  case  where  boundary  layer  growth  is  large,  the 
boundary  layer  thickness  can  no  longer  be  considered  thin  in 
comparison  to  the  radius.   With  these  considerations,  the 
boundary  conditions  for  the  momentum  equation  become 

v  ■>  V(x)  =  A  /(r-6)    as   z  ■*  °° 
This  new  definition  of  V(x)  has  no  effect  upon  the  development 
of  the  u  momentum-integral  equation.   However,  there  is  an 
effect  in  the  v  momentum- integral  equation.   The  development 
is  the  same  as  described  in  Appendix  A  up  to  the  point 


(Id  (vu) 


dx 


|u_  VMi  dr   +  ij^v  dr  \     dz   =   _        /  |v\ 
v  3x        r      dx  r        dx  ;  \  dz} 


Now 


3u        a (Vu)  av 

v  ax  -  dx        -    u  ax 


a  (Vu)    +  u^V    m    /"dr        d6  \ 
ax  r-6        ^dx        dx/ 


Thus    the  momentum-integral   equation  becomes 


f°°  T  a       ,        „    .         2«u    ,      „  •  dr        u-Vdr       u«V  ( &&        dr\*L 
J0  [^    (VU"V^    +  — '<v-V)    E  +  —  &  +  FT  (dx"  "   d^)Jdz 

■■-''BL 
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Simplif ing 

1_  d__ 
r2  dx 


«r2  |"  u(v-V,dz,  +-2r(«.|g) 


u*dz  =  -  v 


(If) 


Introducing  the  velocity  profiles  as  in  Appendix  A,  the  momentum- 
integral  equations  and  the  appropriate  boundary  conditions 
become 


1   d 


c«r  dx 


(r 


•U2-6  [ 

J   n 


f (f-1)  -dn)  +  i!.*|E.U*6 

c  dx 
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i  o 


c*r    dX 
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= 

0 

■^ 

(ii)    • 

vU 
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f" 

= 

u 

—  < 

c 

du 
d"  " 

V2    dr 
c*r   d 

■mt 

(iii) 

f 

= 

1 

1 

(iv) 

f ' 

= 

0 

>  on 

n 

(v) 

f" 
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0 

J 

>  on  n  =  0 


and 


1   d   ,  7  f    .    .   .   "  U«V  1  (&&        6  dr\   f1 

— 2  ""  r2-U-V-6-    f  (l-g)dn  +  — -•- 6    f  dn 

cr^  dx  Jo  r"6  c  ^dX   r  dx/   Jo 


--gg'CO) 

with  the  boundary  conditions 
(i)       g   =   0 


(ii)       g"  =   0 


on  n  =  o 
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(iii) 

g    = 

1 

(iv) 

g"  - 

0 

(v) 

g"    = 

0 

►on   n  =  0 


Again,  it  was  assumed  that  a  basically  linear  relationship 
existed  between  the  boundary-layer  thickness  and  both  the 
momentum  and  mixed -momentum  thicknesses.   Defining  the 
coefficients  Al,  A2,  A3,  and  A4  as  in  Appendix  A,  and  intro- 
ducing a  new  term 

1 

f  dn 
1  -  Dl 


A5  = 


a 


f  (l-g)dn 


D3 


the  momentum-integral   equations    then  became 

I  7VO  ~l  /  \7  \  2 


r '         U1 

§-  +   jj-    (2    +   Al)    + 


2r'         U^        VL 
r  U  V 


6' 
+    6~  + 


A2    ,    r'    /VV 
=   G"+   —'  luj 

1  /V         r'\   R_ 


A4 
G 


The  second  equation  can  be  rewritten 

6«    _   [A4    _   U'    _  V    _   r'    (2    _        A5        \1/Q 
T~       [g~       U~       V"        r"  ^2        r/6    -  lJJ/(1 

Thus,    the   shape    factor  L1    can  be   obtained   from 


+   ,A5    > 
+  r/6    -  1  > 


A2  =  G 


I? 


.  (2  -I-  Al)  + 


(1  -  X?   A3)  + 


L) 


and  the  nev;  values  for  G  and  L^  can  be  calculated  from  the 
value  of  G1,  where  G'  is  now 


(2 11  ♦  V-) 


<r~ 
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and 


G 

L2  ~  2 


["-  -  (tf  f] 


As  the  boundary  layer  grows,  the  axial  velocity  as  defined 
in  APPENDIX  A  for  a  cone  will  no  longer  be  valid.   The  increas- 
ing boundary-layer  thickness  will  cause  the  flow  in  the  axial 
direction  to  accelerate  more  than  would  be  caused  by  the 
convergence  of  the  cone  alone.   The  concept  of  displacement 
thickness  would  take  into  consideration  this  change  in  the 
effective  radius  of  the  cone.   Continuity  now  yields 

Uo'[ro  "  6lx(°)J2,7T  "  U(X)-[ro-(l-X)  -  <51x(X)]2-tt 


thus 


U 


U(X)   =  P 


1  - 


wo) 


-»  2 


(1-X)  - 


*lx(X) 


and 


U 


['♦•^ 


2*  1  + 


(1-X)  - 


6lx(X) 


With  the  new  definition  of  V(X)/  the  conservation  of  angular 
momentum  yields 
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V    •    r      -    6(0)         =      V(X)-    r    •  (1-X)    -    6  (X) 


o        o 


K 


V(X)       = 


r    6  (o ) 

u0-[i  -  t- 


(l-x)    - 


6(X) 
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